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We studied the ultrafast third-order optical nonlinearity in a single crystal of multiferroic bismuth 
ferrite (BiFeOs) in the near-infrared range of 0.5-1.0 eV, where the material is fundamentally trans¬ 
parent, at room temperature. With pump pulses at 1.55 eV, which is off-resonant to the strong 
inter-band charge transfer (CT) transition, we observed instantaneous transient absorption with 
pencil-like temporal profile originating from the two-photon CT transition from the oxygen 2p to 
the iron 3p levels. In contrast, under pumping with 3.10-eV photons, the pencil-like absorption 
change was not observed but decay profiles showed longer time constants. Although the two-photon 
absorption coefficient is estimated to be 1.5 cm/GW, which is ten (hundred) times smaller than that 
of two(one)-dimensional cuprates, it is larger than those of common semiconductors such as ZnSe 
and GaAs at the optical communication wavelength. 

PACS numbers: 42.70.Nq, 75.50.Ee, 78.20.Ci, 79.20.Ws, 78.47.jb 


I. INTRODUCTION 

A lot of efforts have been devoted for seeking materials 
for all optical switching, where large nonlinearity, quick 
response time, and operability at room temperature are 
demanded. Transition metal oxides of Mott Hubbard in¬ 
sulators are known to show ultrafast photoinduced phe¬ 
nomena reflecting broadband optical absorption. Espe¬ 
cially in an one dimensional cuprate of Sr 2 Cu 03 , gigan¬ 
tic photoinduced absorption with an ultrafast decay con¬ 
stant of 1 ps occurs owing to the two-photon absorption 
process [H-Q . Since the wavelength range of probe pulses 
for the process includes the optical communication wave¬ 
length, such kind of materials are hopeful candidates for 
ultrafast all optical switches. So far, one and two dimen¬ 
sional cuprates have been well studied, however, there 
are not so many reports on other transition metal oxides 
especially three dimensional ones. 

Bismuth ferrite (BiFeOs) is one of the multiferroics 
which possesses both ferroelectricity (FE) and antiferro¬ 
magnetism (Tc ^1100 K, Tn ~ 640 K), this material 
is studied from viewpoints of both fundamental physics 
and applications [j, Q. Optical pump-probe measure¬ 
ments under excitation with above bandgap (2.6-2.8 eV) 
photons have been performed in films and single crys¬ 
tals of BiFeOa, where the authors discussed various pho¬ 
toinduced phenomena such as coherent oscillations due 
to magnon and phonon [^-Q, and spectral modulations 
in connection with the transport property [i, [i3. As 
a nonlinear medium, BiFeOa has been studied in terms 
of the second-order nonlinearity in studies such as tera¬ 
hertz emission through optical rectification [HI. and ul¬ 
trafast modulation of spontaneous polarizatio n Il2ll , and 
second harmonic generation spectroscopy [l^ Il4|. How¬ 
ever, there has been almost no study on the third-order 


nonlinearity especially around the optical communication 
wavelength in the near infrared region. 

BiFeOs is also one of the Mott insulators with ultra¬ 
broadband optical response continuously ranging from 
the far-infrared (terahertz) to the ultraviolet regions, re¬ 
flecting strong interactions not only between electrons 
but also between electrons and lattice, spin, orbital de¬ 
grees of freedoms, so that the material makes us expect 
strong and ultrafast nonlinear optical responses in the 
near-infrared below-gap r^ion as ever reported in low¬ 
dimensional cuprates [J, [ 3 . Furthermore, BiFeOs has 
the aspect of a ligand system, in which d-d transitions 
play important role in the optical responses in the visible 
and the near infrared regions, which makes the physics 
of the material more interesting. 

As for the sample preparation, the fabrication of a sin¬ 
gle crystal of BiFeOa has long been difficult. Recently, 
single crystals grown by the flux method have been uti¬ 
lized in many studies [13 [H [Hill , however, the ob¬ 
tained size of the crystals has been limited. In 2011, Ito 
et al. realized the fabrication of large single crystals of 
BiFeOa by using a modified floating zone method with 
laser diodes for heating . Besides the large size (4 mm 
in diameter), the leakage current under high bias voltage 
is extremely low. Hence, we think this crystal will enable 
us to extract only the intrinsic optical property of the ma¬ 
terial by making full use of the interaction length. Thus, 
in the present study, we optically pump the large single 
crystal of BiFeOa with femtosecond pulses and probe the 
transient absorption change mainly in the near infrared 
region, to explore unknown ultrafast nonlinear response 
of BiFeOg. 


2 


II. EXPERIMENTAL 

A large single crystal of BiFeOs with a diameter of 4 
mm was grown with the floating zone method utilizing 
laser diodes as heat sources . By controlling the oxy¬ 
gen pressure during growth, oxygen deficiency is reduced 
a lot so that the leakage current flowing in the crystal 
under high bias-voltage is extremely low. Also the crys¬ 
tal does not show any magnetic hysteresis in the M-H 
curve which is usually observed due to the weak parasitic 
ferromagnetismin in crystals grown by the flux method. 

The unit cell of BiFeOa has the rhombohedral (i?3c) 
structure, which is essentially a cubic elongated along 
the [111] direction. The ferroelectric polarization is along 
this axis. BiFeOa is a G-type antiferromagnet in which 
all the nearest neighbor spin couplings, including both 
inter- and intra-plane ones, are antiferromagnetic. Due 
to the magnetic order, BiFeOa forms the spiral spin struc¬ 
ture with a period of 62 nm along the three equiva¬ 
lent wavevectors of [10-1], [01-1], and [1-10] directions. 
BiFeOa is also a charge transfer (CT) type insulator with 
a bandgap energy of 2.6—2.8 eV, which corresponds to 
the lowest inter-band transition from the oxygen 2p to 
the iron ?>d state. This broad CT band can be further 
decomposed into several peak structures centered at 2.5, 
2.9, 3.2, 4.0, 4.5, and 6.1 eV [U. This CT band with 
strong oscillator strengths brings about strong absorption 
with coefficients of the order of 10® cm“^ [I^. BiFeOa 
also shows relatively weak (^100 cm“^) d-d absorptions 
with broad peaks centered at 1.4 and 1.9 eV, which are 
also split into two by taking the spin degrees of freedom 
into account H- 

Pump-probe measurements were carried out with a ti¬ 
tanium sapphire regenerative amplification system (800 
nm, 35 fs, 1 kHz) as a light source. A (lOO)-oriented 
specimen of BiFeOa with a thickness of 170 fim was 
pumped with the fundamental (800 nm, 1.55 eV) or the 
second harmonic pulses (400 nm, 3.10 eV), and the tran¬ 
sient transmission change was probed with signal or idler 
pulses from an optical parametric amplifier. The dura¬ 
tion of the pump pulse was measeued to be 60 fs, and 
the cross correlation width of the pump and probe pulses 
was 100-110 fs, both at the sample. The penetration 
depths given by the inverse of the absorption coefficients 
[l^ are 33 pm at 1.55 eV, and 0.49 pm at 3.10 eV, re¬ 
spectively. The near-infrared transmission was measured 
with a combination of a monochromator and a cooled 
HgCdTe detector. The 1/e^ beam sizes of the pump 
and probe pulses were 1.35 mm and 130 pm in diam¬ 
eter, respectively. First we pumped the specimen with 
3.1-eV pulses and measured the transient optical con¬ 
ductivity with air-plasma based terahertz time domain 
spectroscopy. As a result, we observed no Drude-like re¬ 
sponse but observed an oscillation of transmission resem¬ 
bling a monocyclic terahertz pulse profile (not shown). 
This result is consistent with the high resistivity con¬ 
firmed with the dc measurements and the oscillation 
profile seems to be caused by the generation of terahertz 


pulses due to the modulation of polarization by ultravi¬ 
olet pulses [H. Since the electrons transport property 
is determined by oxygen deficiencies, the absence of the 
transient Drude response proves the suppression of them. 
Such high quality of the crystal will enable us to extract 
intrinsic property without artifacts. 

III. RESULTS 

Figure 1 shows the time evolutions of the photoinduced 
absorption change at each probe photon energy (0.5-1.0 
eV). The pump photon energy is fixed at 1.55 eV. At 
all probe photon energies, we can see steep increase of 
absorption and its ultrafast recovery. The amount of ab¬ 
sorption change increases with the probe photon energy. 
The inset shows the expanded profiles around the time 
origin. Here the temporal profiles, especially for probe- 
photon energies greater than 0.8 eV, characteristically 
have pencil-like structures with a duration of 100-150 fs. 
This is close to the cross correlation width of the pump 
and probe pulses (100-110 fs). The decay profiles were 
well fitted by the expression assuming two exponential 
functions and a constant as, 

AaL = Ao + Aiexp(-t/Ti) + A 2 exp(-t/T 2 ). (1) 

The constant component (Aq) certainly has a decay pe¬ 
riod beyond 1 ns, however, it did not exceed 1 ms, which 
is the separation period of the laser source. The values 
of T 2 and n are fitted to be 50-60 fs and 1-3 ps, re¬ 
spectively. In the inset of Fig. 2 (a), we present how the 
transient absorption profiles probed at 0.9 eV are decom¬ 
posed into three components. Importantly, the dotted 
curve showing the cross-correlation profiles of the pump 
and the probe pulses coincides with the dominant part of 
the transient absorption change. They appear to be es¬ 
sentially independent of the probe photon energy (wave¬ 
length) within the accuracy of the measurements. In Fig. 
3 (a) we present the dependence of the magnitudes of 
the three components on the pump photon density. All 
the amplitudes of the three components increased with 
the probe-photon energy. More closely, Ai and A 2 were 
proportional to the excitation density, while Aq showed 
the quadratic dependence. We also examined the polar¬ 
ization dependence and found the temporal profile did 
not depend on the polarization of the pump and probe 
pulses, although the magnitude became smaller when the 
polarizations of two pulses were crossed. 

To elucidate the physical mechanism, we changed the 
excitation photon energy to 3.10 eV and measured the 
transient absorption signal at the probe photon energy 
of 0.9 eV. The profile is shown in Fig. 2 (b). The 
decay curve was also decomposed into three parts with 
time constants of 140-170 fs, 1.9-2.6 ps, and nanoseconds 
(which can be regarded as constant). Hereafter, we de¬ 
note the fitting results as Ao(fun) or Ao(sh) to distinguish 
those pumped by the fundamental (1.55 eV) and second 
harmonic (3.10 eV) pulses. The inset shows the profiles 
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FIG. 1: Probe photon energy dependence of the transient 
absorption profiles AaL in a single crystal of BiFeOs. The 
pump photon energy is fixed at 1.55 eV, while the probe pho¬ 
ton energy is varied from 0.5 to 1.0 eV. The excitation density 
is 40 GW/cm^ 



FIG. 2: Pump photon-energy dependence of the photoinduced 
absorption AaL at the probe-photon energy of 0.9 eV. (a) 
Temporal profile pumped at 1.55 eV, and (b) that pumped at 
3.10 eV. Insets show temporal profiles of each (To, Ai, A 2 ) 
component. The excitation density is 56 GW/cm^. Dotted 
curves show the cross-correlation profiles of the pump and the 
probe pulses at each pump photon energy. Here, the profile 
at the pump photon energy of 1.55 eV was directly measured, 
while that at 3.10 eV was estimated as the possible maximum 
one taking into account the group delay dispersion in a BBO 
crystal and optical fiters employed for the measurement. 


of each decomposed decay component. Apparently, the 
decay time of the fastest component for pumping with 
3.10-eV pulses (r 2 (sh)) is longer than that with 1.55-eV 
pulses (r 2 (fun)), where the magnitude of T 2 (sh) is rel¬ 
atively smaller than A 2 (fun). The profile under pump¬ 
ing with 3.10-eV pulses is sharper at the top and does 
not have the pencil-like structure, which appeared under 
pumping with 1.55-eV photons, although the duration 
of 3.10-eV pulses is expected to be longer than that of 
the 1.55-eV ones due to the material dispersion in the 
nonlinear crystal. We also examined the dependence on 
the excitation density as presented with closed markers 
in Fig. 3 (b). A 2 (sh) was quadratically proportional to 
the excitation density up to 60 GW/cm^, approximately, 
while Ao(sh) and Ai(sh) showed linear dependence and 
the slopes become gentle around 40 and 55 GW/cm^, re¬ 
spectively. From these dependences and the absence of 
the pencil-like structure, the origin of photoinduced ab¬ 
sorption under pumping with the second harmonic (3.10 
eV) pulses is different from that with fundamental (1.55 
eV) pulses. 



Excitation density [GW/cm ] 


FIG. 3: Excitation density dependence of the magnitude of 
each decay component under pumping with (a) fundamental 
(1.55 eV) pulses , and (b) second harmonic pulses (3.10 eV) at 
the probe photon-energy of 0.9 eV. For clarity, vertical error 
bars are not indicated, however, the values of error are small 
so as not to effect the dependence. 


IV. DISCUSSIONS 

Now let us elucidate the microscopic origin of the 
present photoinduced phenomena. First, we focus on the 
pencil-like top structure in the transient absorption pro¬ 
file and the origin of the fastest decay component n (fun). 
Ogasawara et al. observed a similar structure in the 
photoinduced transmission change profiles of a quasi-one- 
dimensional cuprate of Sr 2 Gu 03 , which they attributed 
to the two-photon excitation process from the oxygen 2p 
to the even-parity band of Gu which lies at a higher fre¬ 
quency than the one-photon GT band [Il-Q. We plot the 
amplitude of the fastest decay component as a function 
of the probe photon energy as shown in Fig. 4. By over¬ 
laying the absorption spectrum in such a manner that 
the horizontal axis (top) is shifted by the pump photon 
energy of 1.55 eV towards higher energy, the magnitude 
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FIG. 4: Relationship between the magnitude of decay com¬ 
ponents and the linear absorption spectrum of BiFeOs. Aq, 
Ai, and A 2 are represented as closed markers as a function 
of probe photon energy (left and bottom axis). Solid curve 
shows the absorption spectrum (right and top axis) adopted 
from the data shown in Ref. [ 2 J]. 


of the fastest component A 2 (fun) significantly rises up at 
0.7 eV, where the linear absorption with the correspond¬ 
ing photon energy (2.2-2.3 eV) also starts to increase. 
From this and the pencil-like structure, we conclude that 
the two-photon absorption process from the O 2p and the 
Fe 3d level gives the fastest decay component as schemat¬ 
ically shown in Fig. 5 (a). Note the unit cell of BiFeOa 
lacks inversion symmetry at temperatures below 1100 K, 
so that the two-photon transition to the one-photon al¬ 
lowed level is allowed. 

Let us further examine the probe photon energy (wave¬ 
length) dependence. While Ao(fun) significantly in¬ 
creases with the probe photon energy, the increase of 
y4i(fun) is minor. According to the result of the pump- 
probe measurements in the ultraviolet range, the to¬ 
tal lifetime of the CT excited states is in the order of 
nanoseconds M- In the present case, one pump photon 
cannot create the CT excited state, however, the two- 
photon process can realize it. Hence, we consider the 
Ao(fun) component comes from the lifetime of the two- 
photon excited CT states. We can attribute the gradual 
increase of Ao(fun) with the probe photon-energy to the 
dependence of the oscillator strength of the CT transi¬ 
tion on the photon energy. The pump photon energy 
of 1.55 eV is also close to the broadband d-d transition 
from ^Aig to ^Tig, so that the Ai(fun) component with 
the decay constant of 1-3 ps probably corresponds to the 
lifetime of the d-d excited states, and the transition to 
the higher d-d levels or upper Hubbard band induces the 
transient absorption as schematically shown in Fig. 5 
(b). The relatively small amplitude Ai (fun) is due to the 
smaller oscillator strength compared with that of the CT 
transition by nearly one hundred times. As for the chan¬ 
nel of the decay, we can think of some process involving 
the spin system, because the d-d transitions are known 


to be accompanied by magnon excitation in the form of 
sidebands |14|. 


Next let us discuss the photoinduced phenomena un¬ 
der pumping with 3.10-eV photons. The time evolution 
basically resembles the ones reported by Sheu et al. 0, 
who pumped a single crystal of BiFeOa with 3.1-eV pho¬ 
tons and measured the transient reflection change also 
at 3.1 eV. From the quadratic dependence of A 2 (sh) 
on the excitation density, some cooperative interaction 
is certainly involved in the decay process. Similarly to 
the case in the previous report, the relaxation from the 
CT excited states to the bottom of the band through 
electron-electron or electron-phonon scatterings is highly 
possible as the origin. As already written, T 2 (sh) slightly 
decreases with the excitation density, which seems to sup¬ 
port this speculation. Furthermore, while Ao(fun) shows 
the quadratic dependence up to the excitation density of 
56 GW/cm^, Aq (sh) shows the linear dependence only 
up to 40 GW/cm^ and then the slope becomes gentle. In 
the case of excitation with fundamental (1.55 eV) pulses, 
some cooperative interaction seemingly works between 
the two-photon excited CT states. This is the origin of 
parabolic dependence of Aq on the excitation density. 
In contrast, Ao(sh) reflects the one-photon excited CT 
states, it can easily be saturated at the relatively weak 
excitation density as we observed in Fig. EKb). 


Now let us discuss the two-photon absorption coeffi¬ 
cient of /3 in BiFeOs. Using the transient transmission 
change of -0.11 at the probe photon energy of 1.0 eV 
(1.24 /rm), P is estimated to be 1.5 cm/GW. The P in 
the one-dimensional cuprate of Sr 2 Gu(^ is 160 cm/GW 
at the peak of the two-photon bandQ, while that of 
the two-dimensional Sr2Gu02Gl2 is 12 cm/GW. This di¬ 
mensionality dependence of the optical nonlinearity in 
materials relating to high Tq cuprate superconductors 
was closely examined by Ashida et ai, in which the p 
in the two-dimensional cuprates was typically ten times 
smaller than that of the one-dimensional ones. The mech¬ 
anism was explained using the cluster modelQ. Intu¬ 
itively, to limit the coherent vibrational oscillations of 
electrons with large amplitudes in one direction can en¬ 
hance the optical nonlinearity. While in the case of two or 
three dimensional materials, motion of electrons spreads 
into multiple directions even if electrons and lattice are 
vibrated in one direction by linearly polarized optical 
pulses. Thus the low-dimensional cuprates may be more 
beneficial in terms of application to optical device. How¬ 
ever, the p value is still larger than those of common 
semiconductors such as ZnSe and GaAs at the optical 
communication wavelength [2^. BiFeOa also has the 
strong aspect as a multiferroic, multifunctional responses 
relating to the magneto-optic effect can be expected in 
low temperature or polarization sensitive measurements. 
To explore such phenomena is our future work. 


Finally, let us give a brief comment on the relationship 
between the present results and the optical anisotropy. 
BiFeOa shows uniaxial birefringence with an extraordi¬ 
nary axis along the [111] direction. However, we observed 
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(a) (b) 



FIG. 5: Schematic drawings of the electronic states of bismuth 
ferrite to show the origin of the ultrafast decay components 
under pumping with fundamental (1.55 eV) pulses, (a) Upper 
Hubbard band of Fe 3d and O 2p level, in connection with the 
fastest (H 2 (fun)) and the slowest ylo(fun) decay components, 
(b) Relationship between the d-d levels and the second fastest 
(Hi(fun)) component. 

no significant polarization dependence in the results of 
the pump probe spectroscopy. We can understand the 
reason considering the three dimensional chemical bond¬ 
ing structure in BiFeOs; in the present phenomenon, the 
CT excitation and d-d transition play important roles, 
while absorption due to phonons and magnons shows 
significant polarization dependence. At low tempera¬ 
tures, the coherence of elementary excitation will grow 
up enough to defeat dephasing factors and we will ob¬ 
serve the transient anisotropic optical response originate 
from coherent phonon and magnon. 


of BiFeOs using sub-lOO-fs optical pulses in the near- 
infrared and ultraviolet ranges. With pump pulses at 
1.55 eV, we observed sharp decrease in transmission 
which remained constant around the minimum for ap¬ 
proximately 100 fs, it decayed with a time constant as 
short as 50-60 fs. Examining the dependences of the 
signal intensity on pump fiuence and probe photon en¬ 
ergy, we conclude the dominant transient absorption orig¬ 
inates from the two-photon absorption process from the 
oxygen 2p to the iron 3c? level. The decay profile also 
has two other components with constants of 1-3 ps, and 
nanoseconds which can be fitted as a constant in the 
present observation range. We attribute them as some 
relaxation process involving spins, and the total lifetime 
of the CT excited electrons via the two-photon process. 
The two-photon absorption coefficient is estimated to be 
1.5 cm/GW, which is ten (hundred) times smaller than 
that of two(one)-dimensional cuprates. This result shows 
such ultrafast response is common in strongly correlated 
systems and may reflect the dimension dependence of the 
nonlinearity. The present /3 value is still larger than those 
of common semiconductors such as ZnSe and GaAs at the 
optical communication wavelength. With pump pulses at 
3.10 eV, the pencil like two-photon absorption did not ap¬ 
pear, but the profiles had longer time components of 140- 
170 fs, 1.9-2.6 ps, and nanoseconds. We attribute them 
to the collective electron-electron and electron-phonon 
scattering in the GT band, some decay process involving 
spin or band to band relaxation, and the total lifetime of 
the GT excited levels, respectively. 


V. CONCLUSIONS 

In conclusion we demonstrated the ultrafast pho- 
toinduced absorption spectroscopy in a single crystal 
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